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ABSTRACT 



We derive the closed form expression for the bit error probability of asynchronous 
dense WDM systems employing an external OOK modulator. Our model is based upon a 
close approximation of the optical Fabry-Perot filter in the receiver as a single-pole RC 
filter for signals that are bandlimited to a frequency band approximately equal to one 
sixtieth of the Fabry-Perot filter's free spectral range. Our model can handle bit rates up to 
2.5 Gb/s for a free spectral range of 3800 GHz and up to 5 Gb/s when the power penalty 
is 1 dB or less. 
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I. INTRODUCTION 



Asynchronous wavelength division multiplexing (WDM) systems have been 
increasingly proposed as an attractive alternative to coherent optical frequency division 
multiplexing (FDM) systems [Ref. 1-5], Although asynchronous WDM systems with 
direct detection do not have the channel capacity of coherent optical FDM systems, they 
are much less costly to implement. Furthermore, present filter technology enables the 
designers to tightly pack the channel, resulting in asynchronous dense WDM systems that 
can provide aggregate bit rates of many terabits per second (1 Tb/s = 10'^ b/s). 
Asynchronous dense WDM systems are particularly attractive in the area of undersea 
surveillance where hundreds of sensors and data collection sites are envisioned being 
merged onto single-fiber superhighways through massive data fusion. Other applications 
call for relatively low-bandwidth data collection over many months to be dumped quickly 
to a remote recording site in a matter of minutes. This "collection-and-dumped" 
compression can demand total data rates on the order of hundreds of Gb/s. Long distance 
between data collection sites and a remote recording site requires the use of optical 
amplifiers. Therefore, it is necessary to pack all channels within the optical amplifier 
bandwidth. 

An asynchronous dense WDM receiver with on-off-keying (OOK) modulation can 
be modeled as shown in Fig. 1 . Conceptually, the analysis involves two main operations: 
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1) a convolution operation to evaluate the signal at the output of the optical filter, a 
Fabry-Perot (FP) filter in our investigation, and 2) the integration of the output of the 
photodetector. Evaluation of bit error probability by the numerical analysis of these two 
operations has been carried out in [Ref 6], with a number of approximations made to 
reduce the computational complexity. In this investigation the FP filter is shown to be 
well approximated by an RC filter within the frequency range \f-fg\< FSR / 20tc, where 
FSR is the free spectral range of the FP filter [Ref 7] and /q is the FP filter center 
frequency. For example, given FSR = 3800 GHz, the approximation works very well for 
I f-fo I < 60.5 GHz; that is, the effects of adjacent channels within a 121 GHz bandwidth 
centered at fg must be included, while all others can be neglected. This simple model 
agrees well with [Ref 6] as demonstrated in Section HI. Furthermore, this model enables 
us to obtain a closed form analytical expression for the bit error probability for which 
numerical results can be obtained with little effort. Our investigation shows that this 
simple model provides accurate results as compared to those in [Ref 6] for bit rates up to 
2.5 Gb/s when the effects of two adjacent channels are included with FSR = 3800 GHz. 
Actually, when the power penalty relative to single channel operation is 1 dB or less, there 
is virtually no difference in the effect of four or two adjacent channels. Thus, for this 
power penalty criterion, this simple model can handle bit rates up to 5 Gb/s for a FP filter's 
FSR = 3800 GHz. 

In Section II the closed form expression for the decision variable, and consequently, 
the bit error probability assuming all channels are bit asynchronous as in [Ref 6] is 
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derived. Section III presents the numerical results which include the bit error probability 
versus the signal-to-noise ratio as function of the FP filter bandwidth and channel spacing, 
and the power penalty (relative to single channel operation without filtering or with 
filtering but no intersymbol interference) versus the channel spacing as a function of the 
bandwidth. Finally, a summary of results appears in Section IV. 
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n. ANALYSIS 



The receiver model for the asynchronous dense WDM system is shown in Fig. 1. 
The desired signal is filtered by a Fabry-Perot (FP) filter that rejects adjacent channels. 
The photodetector is assumed to have a responsivity St (AAV). The detected current is 
amplified by a low noise amplifier that contributes a postdetection thermal noise n(t) with 
spectral density Ng (A^/Hz). The decision variable at the output of the integration is 
compared to a threshold a to determine whether a bit zero or bit one was present. 

A. INPUT SIGNAL 

For convenience, we designate channel 0 as the desired channel, and channel k as 
an adjacent channel where k = -M2, ..., -1, 1, ..., M/2 with M an even integer. We 
consider the equivalent lowpass (complex envelope) data signal in channel 0 and channel k 
as follows: 

^o(0 = Z bQjPT{t-iT) (1) 

/=-Io 

bk(t) = E ® 

l=-L 

where 

T: bit duration 

{0,1 }; bit in channel 0 in the time interval (/T, (/+1)7) 
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e { 0, } is the /"■ bit in channel k in the time interval {IT + Xjt, (/ + 1)7’+ Xk) 

(j)^ : a phase offset between channel k and channel 0 and is assumed to be uniformly 

distributed in (0, 27t) radians 

©k : radian frequency spacing between channel k and channel 0 with = -©^ 

Xk ; is the time delay between channel k and channel 0 and is assumed to be 
uniformly distributed in (0,T). 



The function P^t - iT) is defined as 



Pr(?-<T)={i; 



/7’</<(/+l)7’ 

otherwise 



( 3 ) 



In both (1) and (2), the non-negative integers and L represent the number of bits in 
channel 0 and k, respectively, that proceed the detected bits b^Q. The received 
asynchronous dense WDM equivalent lowpass signal at the input of the FP filter is given 
by 

MU r— 

r{i) = ^[Pboit)+ X yfPbk(t) (4) 

k=-M/2 
k^O 

where P is the received optical power. 
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B. FABRY-PEROT FILTERED OUTPUT SIGNAL 



The FP filter can be characterized by the following equivalent lowpass transfer 

function [Ref 1,7] 

H( = ^~P ^ l-A-p 

) j_pg-;2jy/KSK i_p 



1-p 



l-pcos(^H/psin(^) 



1-^4 -p 

1-p 



(5) 



where p is the power reflectivity, A is the power absorption loss (zero for an ideal FP 



filter) and FSR is the free spectral range. For \ f\< FSRJIQti and assuming .^ = 0, we can 
approximate // (/) as follow: 



^(/)- 



1-p 






1 



1 +/ 



(1-p)/^ 



1 




\f\<FsmoTi 



(6a) 



where 



F5R(l-p) 

P 



(6b) 



The free spectral range FSR can be related to the full width at half maximum (FWHM) 
bandwidth B and the finesse F of the FP filter as 



FSR = 



1-p 



= BF 



(7) 
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Thus if the signal is bandlimited to \ f\ < FSR/20n, we can truly approximate (5) with a 
single-pole RC filter with the following transfer function and impulse response 

H{f) = 7 ^, ( 8 ) 

1 +/— 

h{t) = ce~^‘ , t>0 (9) 

Figures 2a-b show the magnitude and phase (radians) o^ H{f) of the FP filter in (5) and 
its single-pole RC filter approximation given in (8) for p = 0.99, F = 312.6, B = 12.16 
GHz and FSR = 3800 GHz. Note that as the frequency increases, the phases of the FP 
filter and the RC filter differ markedly, but the magnitudes of their transfer functions 
remain identical and attenuate rapidly. When | / | > FSR/IOk, the magnitude of H{f) is 
very small, and therefore, the effect of adjacent channel interference beyond this frequency 
range is negligible. Fig. 3 shows the normalized impulse response of both FP and 
single-pole RC filters. In summary, the above approximation is valid for asynchronous 
dense WDM analysis when the filter finess F is large or equivalently the FWHM 
bandwidth B is small since the equivalent lowpass signal must be bandlimited to about 
I f\<FSRI2Qn. 

This approximation has been used in [Ref 5] to study spectral efficiency of optical 
FDM/ASK systems, which involves the evaluation of the decision variable for worst-case 
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analysis using the eye diagram technique. Since we are interested in the detected bit b 
in the time interval (0,7'), we consider the output filtered signal 5(/), Q<t<T given by 

5(0 = SB{t) + Sisi{t) + SACl{t) , 0<t<T (10) 



where 

Sg (/): desired signal 

Sjgjit ): intersymbol interference (ISI) signal 
Sac, (0- adjacent channel interference (ACI) signal 
These signals are evaluated using (4) and (9) as follows (detailed derivations in Appendix 
A): 

t 

SB{t) = JPbo^o 1 h(t-X)cDi 

0 

= 0<t<T (11) 




= Z (12) 

;=-Io 
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SACliO 






r— MU 

= JP 2 

k=-MI2 

M 


_2 (/+I)7-+Xi 

[2 At,, J 

’ /r+Tifc 


o<f<r 




+ At,-i j A(/ - 

-r+Tfc 


0 < ? < T* 




+ At,-1 J 

-T’+Xx, 


%k<t<T 




+ At,o 1 A(< - ] 

Xifc 


Xk<t <T 



= JpT [i 7t(0, T)bki e-" g(c+;<o»)(;7'+x»)(g(c+y<o,)7' _ 

k=-MH l=-L ’ l+j'-T 

k*0 

+n(0, Xk)bk,-t 7 ^ - g(c+;cofc)(-r+Tfc)^ 

l+/~ 

+7t(Tt, 7)6*:-i e-" 

i+;-r 

+7t(Tt, 7)6i,o e“" ^ ( 13 ) 

1 +/ c 

0<t<T 



where 






l,/l</</2 
0, otherwise 



is a unit pulse between t, and 
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The FP filtered output s(i) is detected by the photodetector which produces a 
current of d \ s(t) p Amps. This current plus additive white postdetection thermal noise 
current from the amplifier is integrated by the integrator to obtain a decision variable for 
the threshold detector. 

C. DECISION VARIABLES 

The decision variable V appearing at the integrator output consists of the signal 
component X and noise component N 

Y = X+N (14) 

where 

T 

X=\ As{t)\'^dt (15) 

0 

T 

N =\ n{t)dt (16) 

0 

We note that V is a zero mean Gaussian random variable with variance N^T. Substituting 

(10)-(13) into (15) we obtain the signal component Xas a function of the three parameters 
Xk, cT, and which represent the effect of intersymbol interference and adjacent 
channel interference. 

i=~-Lo 
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M/2 -2 



e~J'^k'’k 



2cT 



k=-M/2 l=-L 

k*0 



i+y? 



cT 



M /2 

I 

k=-M/2 

k*0 


M /2 

S 

m--MI2 


bk.-xbl^., e-J^'-k-k-^m^m) 


rCTj^; COJt— CO/n -I 


(i+/?-)(i-;^) 


/ \ / CDj^T^CO ffi 


. M /2 

: S 

k=-M!2 

k^O 


M /2 

I 

nt=-MI2 

nv^O 


e-J^k^k-^n,Xm) 




(i+;^Xi-;^) 


[ l+;i^ 



g(c-j^mX-T+^m) (e~^^^—l) 

h 5 ^ ' 



^ pj M/2 M2 

+ ^ jt; 

k=-M!2 m=-M2 2 (l+ 7 -^)(l-;^^) 

fcto m^O 



(c+7<0 *)TA: ^ 2 _ ^ -(c+;(0 it)7’ (c-;CO m)Tm ^ J ^ -(c-> >0] 



_|_ ^ P_l_ ^ ^ ^l-CZk (Ofc=i 



CO, 



cT 






k=-M/2 m=-M/2 {^+J-tA^-J-r) - — — ;^ (Hk^a^m 

k*0 m*0 J(<s>^r^m)lc '” 



^ PT bkfib*„Q e ^ g^^j'^0'^k gHc^<^>m)-ig'^ 

^ ~ ^ u7 TT [ 



k=-M/2 m=-M/2 (1+7^)(1-;^) 

JTT^O 



y+j^ 



gic-ja>m)tm^g-<o-J<^k^T_ ^-(c-jo)k)xg'^ (g-2c7’_g-2ctg^ g(c+;tojt)t|. g(c-j<i>m)im ^ 

h¥ ^ ^ 
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•7 PT 

+ ie^Re{ Z 

fotO 



m=-M/2 

m^O 






L I . .tdm J / 



I+y^ 



+ ^^Re{ S S 2; ^ 

k=-M/2 m=-MI2l=-L (i+j-rXH^) 



[ 2 



-.pT M/2 M/2 -2 

+ /f 2 S S 

yt=-M2 m=-M/2 l=-L 

h^O m^O 



g-j(o>kVlc-o>mXm)(^g(c+ja>ic)T_ ^(c+ja>kXlT+^0 

(i+;?-Xi-y™) 



[ 






+ 






1+7^ 



2 



]} 
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MU 

Z 



MU 

Z 



^(ci) X m “C m ) 



k=-MH m=-MH ( 1 +/“)( I ~j~T') 

fctO m*0 



[ 






1-y^ 



2 

bic-ib^, 

_ (cat-um) 

; — s — 



6a: _1 A * 0 







l+j^ 



_l j I 



(The above whole term will be zero if © = © ;« ). 



i=-Lo 
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+ /P—lfKe { 2 ^ 2 ^ — ^ ^ -. 

/ t =- M /2 /=-L 2 ( l + j ^) 

k^O 



W bofibk,-i e-i'^k^k e>'^k^k-\ 

' " , .“/t L “;t 

k=-MH 1+y-r J-r 

fcto 



“wf 2 



] 



^ bofibk-\ e~J'^k^k 

■*■ ^ 7~~^ L 9 J 

k=-Ml2 1+y-r ^ 



^2 bofibk.Q e ^“‘k'^k gj<^k'^_gj<^k'^k 

" l “ 2 ^ ml L ' 



0)1, 

k=-M/2 l+J— 

kio 












]} 






-1 



cT ,. 



/■=-Z,o 



Re { L L — ST 

k=-M!2 l=-L 2(1+;-^) 

k^O 



M2 e-j<^kH-e-^k 

" I " 2 < (oi. L ' 



0 ) 7 , I 

k=-M/2 \+j-r ^-j-T- 

fatO 



( e -^ k - e ^ ke -^^ T ) ( i - g -^^ J '^ k )^) 



] 



M2 a. 

+ I " 



0 re 



~‘^'^k-e-J“‘k^ke-<^J<‘k)'^ e~‘^'-k-e‘^^k g-2c7' 



k=-M/2 1+y? 

/ bK ) 






]} 



(17) 



where = min (xjt, Tm) "rg = niax(T/c, Tm) 
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D. BIT ERROR PROBABILITY 



For a detection threshold a and an ISI/ACI bit pattern b = {b^ , b^^ 

/ = -L, 0; k= -M/2, M/2, 0; the conditional bit error probability of the OOK 

signal represented by the Gaussian random variable Y in (14)-(16) is given by [Ref. 8] 



Pe(b) = 7S(^w) + \Q (^w) 



a-Xoib), 



Mr 



Mr 



where Q (X) is defined as 



(18) 



Q(x) = I e 



(19) 



and Xq and X^ are the values of A' in (17) for = 0 and = 1, respectively. The 
average bit error probability is obtained by taking the expected value of (b) in (18) 
over all bit patterns b. The minimum bit error probability is obtained by optimizing over 
the threshold a. In summary, given a, we calculate the following expectations: 



Pe =E {Pern 

b 



( 20 ) 



P e, min =rmii£’ {Pe(b)} 



2M(L+\)+Lq 



I p(b) 

2^ patterns 



15 



Pib ) = f ( ^ 



In T 

[J-J J-Je( 

MO MO 



a-Xo(^^i2 . . .4>m/2 




-MU • • • d^MH dx-Mi2 • • • d%ua 



2n 
M 0 



■Jec 



hi 0 



X 1 ^M/2 )~Q 

JW 



)d(^-M/2 • • -<^<))A//2 dX-M/2 • • .<^XA//2] 



( 21 ) 
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ra. NUMERICAL RESULTS 



In this section we present the numerical results for a) bit error probability versus 
signal-to-noise ratio Z = ^P^TINo as a function of the normalized channel spacing 
(normalized to the bit rate) 



I = ^ = (Af)T ( 22 ) 

where A/'=co^/27tA is the equal channel spacing in Hz, and b) power penalty versus 
normalized channel spacing / as a function of the FP filter parameter cT={%l Jp)BT 
where B is the filter full width at half maximum bandwidth(FWHM). For our model to be 

valid we set M = 2, that is, we constraint the signal to be bandlimited to within two 

adjacent channels, thus we incorporate the degradation the two nearest adjacent 

channels. We observe that for bit rates of 1 Gb/s or less, oui model is valid up to A/ = 10 

and very little difference is observed between M= 2 and M= 10. Also, we observe that 
there is little difference between M=2 and A/= 4 when I>10 for bit rates up to 3 Gb/s. 
In all results we set Lg = 2 and L = 1. 

Figures 4-5 show the minimum bit error probability versus signal-to-noise ratio 

( Z ) for c7’= 5 and 10, respectively. Along with each curve, we also show that of the 
synchronous case, i,e, Xa = 0, = 2 and L = 0, and that of a single channel (SC) 

operation without filtering or with filtering but without ISl. In Fig. 4 we observe that a 

large degradation occurs due to ISI for cT= 5 which represents a narrowband filter. As 
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the FP filter bandwidth is made larger as in Fig. 5 with cT= 10, the ISI is reduced but the 
ACI increases. 

In our model, we are constrained to M = 2 for the case under consideration. We 

use an FP filter with FWHM 5=12.16 GHz, free spectral range FSR = 3800 GHz, finess 
F=7cyp /(I -p) = FSR/B = 312.6, and c = 38.4 GHz; 1/T = 2.56 Gb/s, and cT = 15. If 
the power penalty related to single channel operation is to be less than 1 dB, then by 

Fig 6, a minimum normalized channel spacing of / = 12 must be used. Equivalently 

from [Eq. (22)], we could use 6f = IfT = 12*2.56 GHz = 30.4 GHz, (i.e. the channel 

spacing is a multiple of the bit rate). In this model the farthest adajcent channel for M= 2 

is twice the channel spacing which is 60.8 GHz. This verifies the assumption I f-f^ I < 

FSRUQtz = 60.5 GHz, where fg is the FP filter center frequency. This result agrees well 

with that in [Ref. 6; Figs. 6,9, M/F = 0.4, a = 0.2]. Thus we incorporate the degradation 

caused by the two nearest adjacent channels. We observe that for bit rates of 1 Gb/s or 

less, our model is valid up to M = 10, and very little difference is observed between M = 

2 and M = 10. Also, we observe that there is little difference between M = 2 and M = 4 
when / > 10 for bit rates up to 3 Gb/s. In all results we set Lg = 2 and L = 1. 

Figure 6 also shows the power penalty for an asynchronous dense WDM system 

relative to a single channel operation at the minimum bit error probability of 10''^. This is 

the required additional signal power (dBW) for the asynchronous dense WDM system to 

be able to operate at the 10''^ bit error probability achieved in the single channel system 

with a SNR=l2dB. The asynchronous dense WDM system is /57-limited at 2. 15 dB, 0.95 

dB, 0.6 dB, and 0.5 dB in power penalty for cT = 5, 10, 15, and 20, respeaively. It is 
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seen that for a 2.3 dB power penalty, the normalized channel spacing can be as close as I 

= 6 (i.e., a channel spacing of six times the bit rate) for cT = 5. If the power penalty 
criterion is 1 dB, the normalized channel spacing is/ » 12 for cT= 10, 15, 20. We remark 
that although the exact transfer functions of the FP filter is used in plef 6], a number of 

approximation have been made to obtain numerical results. The approximations are 1) 

the ISI is obtained by modeling FP filter as a single-pole RC filter [Ref 6, Eqs. (4) and 

(36)], 2) approximating the finite integration with an infinite integration in the calculation 

of ACI [Ref 6, Eq. (15)], and 3) the beat interference is ignored. On the other hand, the 

ISI and ACI in our investigation are obtained by modeling the FP filter as a single-pole RC 

filter, using finite integration and including the beat interference. Since the results in our 

investigation and in [Ref 6] agree well, we conclude that approximations are quite valid. 

We also note that our results also agree well with the simulation carried out in [Ref 1, 

Fig. 17], 

The above numerical results shown in Figs. 4-6 are obtained with an optimized 
threshold setting. Figure 7 shows the power penalty for fixed threshold a = £PT!2 which 
is the same optimum threshold for single channel operation (midpoint between the 
received power for bit zero and bit one). It is seen that the performance of an 
asynchronous dense WDM system is quite sensitive to a for a narrow band filter. An 
additional 1.5 dB is observed for cT = 5 for / > 8, and 0.4 dB for cT = 10 for / > 12. 
Negligible degradation is observed for cT= 15, 20 for/> 16. 
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Figures 8-9 show the power penalty versus normalized channel spacing as a function 
of FP filter parameter cT for the worst-case analysis with optimal threshold and fixed 
threshold, respectively. The worst-case bit pattern is fixed to produce the minimum X, 
and maximum A'o where X, andXg are the values of X in Appendix- A equation (4-21) 
with ^ 0.0 “ ^ ^ 0.0 “ 0 respectively. 

We observe that the power penalty for the worst-case analysis is only slightly larger 
than that of the exact analysis for /> 10 shown in Fig. 6. Similarly the power penalty for 
the worst-case analysis with fixed threshold is only slightly larger than that of the exact 
analysis with fixed threshold for / > 10 shown is Fig. 7. The reason for this is that for 
large channel spacing (/ > 10), the ACI effect is small, so the ACI bit pattern has a small 
influence on the power penalty. 

Figure 10 shows the normalized optimal threshold for the exact analysis shown in 
Fig. 6. It is observed that a »0.4 for / > 10. Note that the normalized optimal threshold 
for the single channel operation is a = 0.5. 
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rv. CONCLUSIONS 



We have presented a simple model for the analysis of asynchronous dense WDM 

systems employing an external OOK modulator. The only approximation that we use 

involves the modeling of the Fabry-Perot filter by a single-pole RC filter assuming the 
equivalent lowpass signal is bandlimited to the frequency range I/I <FSRIlQii. This 
model enables us to obtain a closed form expression for the bit error probability which 

previously could only be obtained via numerical analysis [Ref 6], For an FP filter with an 

FSR around 3800 GHz, our model can include the ACI effects of two adjacent channels 

for bit rates up to 2.5 Gb/s. Our numerical results show that this model agrees well with 

that in [Ref 6]. 
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POSTDETECTION NOISE 
n (t) 



Figure l:OOK receiver structure 
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Frequency (GHz) 

(b) 



Figure 2; Spectral characteristics of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 
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Figure 3; Normalized rimpulse response of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 
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Z(dB) 



Figure 4; Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for c7=5 
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Probability of Bit Error 




Figure 5: Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for cT=\0 
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Power Penalty (dB) 




Figure 6; Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT 
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Power Penalty (dB) 




Figure 7; Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT with a fixed threshold a=0.5 
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Figure 8; Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT for worst-case with optimal threshold a={XQ+X,)!2 
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Figure 9: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT for worst-case with fixed threshold a=0.5 
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Optimal Threshold 




Figure 10: The normalized optimal threshold versus normalized channel spacing as a 
function of Fabry-Perot filter parameter cT 
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APPENDIX A 



Derivation of Formula 



L INTRODUCTION 



Desired channel Channel 0. 

Adjacent channel : Channel k. 

k = - M/2, ... -1, 1, ... M/2, M : even. 

Bit in channel 0 in time interval {iT, (i+l)T) : b^ . e { 0, 1 } 

Detected bit in channel 0 in (0,T) : 6 { 0, 1 } 

Bit in channel k in time interval (IT + Tk, (/ + 1)^ + Tjt) ' 
bk,i e 

frequency spacing between channel k and channel 0, ©a: = — cO-A:. 

(|)^ ; phase offset between channel k and channel 0, uniformly 
distributed between (0, 2 tc). 

Xk ' time delay between channel k and channel 0, 0 <t*<T. 



0 

Data signal in channel 0 : 6o(0 = S ^o,/ PT(t — iT) 

i=-Lo 

0 

Data signal in channel k : ^ bkje^^~'^'‘^ Pr(t- (IT + Xk)) 

l=-L 

Lo, L : positive integers. 

P (i\ = / 

T\}) X O^otherwise 

The received signal at the input of the FP filter of channel 0 is : 

_ MU 

r(t)=JPba(t)+ E 

k=-MH 

fcto 

P : received optical power. 




( 1 - 1 ) 

( 1 - 2 ) 
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The output of the FP filter is: 



00 



ro (0 = J h{t-X)r(X)dX 



—00 



1 — 1 ' /— 7 

= Jp\ hQ-X)bo(X)dk+JP £ J hil-X)bk(X)dk 



—00 



k=‘M/2 -00 



00 



= -IPbofij Kt-X)PT(X)dX 



—00 

-1 “ 



+ Jp't bo,i j Kl - X)Pt{X - iT)dk 



/— Z/Q —00 
MU -2 



00 



mi A -z c X 



k=-MH l=-L -00 



00 



+i>t -1 j AO - - (-r+ 

-00 

00 



( 1 - 



-00 



is the equivalent lowpass impulse response of the FP filter of channel 0. 
Since the detection interval is 0<t<T, we only need to evaluate 

S(t) = (t), 0<t<T and write S(t) as : 

S(t) = S,(t) + SJt) + S^ci(t) o<t<T 

Sg (t) : Desired signal in channel 0, the 1st term (i=0) of (1-3) 
Sjgj(t) : Inter-symbol interference, the 2nd term (i ^-1) of (1-3) 
Adjacent channel interference, the 3rd term (Z) of (1-3) 

k 

All evaluated in 0 <t<T. 



33 



II, CALCULATE (t) and S ,., ft) 



In this section we are to calculate the easy part of S(t) — Sg (t) and Sj^j (t). 

A. FP filter ; 

The lowpass equivalent transfer function of FP filter is 

1-p 1-^-p _ 1-p l-A-p 

^(f)^ 1 _np-j2n/IFSR • 1-0 ~ .2nf. . . ,2nf. * l_n 

^ ^ P l-pcos(^)+ypsm(^) P 

p : power reflectivity. 

A power absorption loss ( 0 for ideal filter). 

FSR : fi'ee spectral range. 

Since /« fsr (for operating frequency range), we can approximate 
H (f) 2 iS (assume A = 0) : 



m = 



1-p 

(i-p)+/p 



2nf 

FSR 



!+/•■ 



1-p 

27t/p 
(l-p)F5R 



1 

1 -271/ 

1+/-T 



where c = 



FSR{\-p) 



for FP filter, we also have : 



ron 7C Jp 

, B : full width at half maximum bandwidth (FWHM) 
of half power bandwidth. 

Using inverse Fourier transform, we find 



h(t)= { 



ce t >0 
0, otherwise 
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B. 



Calculate (t) : 

00 

S^(t) = Jp Z?o,o J h{t — X) Pt{X) dX 



-00 



i 

= JPbo^o I h(^ - X) dX 



i 

= -JPbofi J dX 



= JPb 



t 



0,0 ce j dX 



0 



= JPbo,o(\- ^-"0 



(b) For Bit "0" : (0 = 0 



C. Calculate Sj ^j (t) 



-1 j. 

Sisi (0= Jp S ^0,; f h(t - X) Pt (X - iT) dX 

i=-Lo 



00 



-1 



-00 

(;+l)r 



= jp H bof f h{t-X)dX 

i=-U 



= Jp 2 ^ 0 / f ce ^0 

/=-Z,o 



iT 

0 + 1 ) 7 ’ 



iT 



-1 



i=—Lo 



0<t<T 



(a) For Bit "1" : S^j (t) = JP {\ - 0<t<T 



0<t<T 



= JPe~^^ 2 bo,i(eO+^'^^^-e^^^)dX 0<t<T 



( 2 - 



( 2-2 
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II L CALCULATE S^ ft) 

In this section we are to find with the procedure : 

A. When/<-2 

B. When / = - 1 

C. When / = 0 

D. Summary for (t) 

A. When / < -2 

Integrate the first term of b^(t) in (1-3), we have : 

M2 -2 “ 

Jp I. 2 btj f h(t - X) - (IT+ Ti))^a 

k=-MH l=-L _io 

_2 (/+1)7’+Xfc 

= JPZtbk,i J h(t - X) 

^ ^ IT+Xk 

_2 (M)T+Xk 

= JP HZ bk I J 

" ' ’ /r'x, 

_2 (/+i)r+Tit 

= ^fP HZ bkj J 

^ ^ /r+xfc 

= JP22^ Q fcXfc)g (c+;co fciC/r+Xfc) (c+ 7 'co jfc)r _ | ^ 
k I 

= JP 22 bicj e“" ^ e(^+/“»X'7'+T,)(g(c+>t)7- _ J) 
k I ^+J— 

= JP 2j 2 Q-<^{t-'^k)^Q(c+j(Sik){lT+T) _ g(c+/cofc)/7’^ 0</<7’ 

A: / 1+/— 



(3-1) 

(3-la) 



36 



B. when / = -1 



B-1 / = -l, 0 < / < r 



Integrate the second term of (t) in (1-3) : 
MU 



00 



■JP S bk.-\ I h{t - X) -(-T+ xt))dk 



k=-MH _oo 



I 

= /PHbk.t I h(t-X)eJ“'‘^^-''=^dX 



-T+Xk 

t 



= JP'Z bk-i I 

-r+tt 



= 7? S *t,-i e-" 






^ i"y c 



0<KTjt 



g-2 / = -l, Xk<t<T 

MU . 

JP E 6t.-i J /i(/-A,)e>‘(^-''‘>PK>‘-(-7’+xt))ca 



00 



k=-M/2 

htO 



-00 



Tit 



= JPZbk-i I /!(<->.) 

^ -r+Tfc 

e~}'^k'’k 

Dk-\ e - 
k 

b. 



= Jp'Lbk-\ e-“ ^^£<‘^■'“‘>’*(1 

l+y-r 



= 7 ? E ( 1 - £-<«>*)?■) 

^ i"y c 



Tk <t <T 

Tk< t <T 



(3 

(3- 
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C when 1 = 0 

C-J. I = 0, 0<t<Zk 

CX3 

■JP S bkj, J h(t - Pt-(^ - Xi) = 0 0 < / < Xt ( 3 - 4 ) 

/t ’ jL 



C-2 / = ^, Tk<t<T 

Integrate the third term of (t) in (1-3) : 



JP'Lbkfi J 

k ’ ^ 

A. — oo 

t 

= JPl.bko f *(/- dk 

* x» 

= Vp Z e-"' ^ 

1+y-r Xk<t<T 

= Jp ^ — M- _g-C(r-TAr)^ Xk <t <T 



(3-5) 

(3-5a) 



A Summary for 

Define : Tl{t^ , = f/,y (t> - f/,2 (t), t>0, t\<ti 

where UJt) = {\^^ c>0 

is the unit step function. 
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then 71 (tj, ( 2 ) — otherwise 

is a unit pulse between tj and t 2 . 

hence for 0 <t <T (t) can be written as ; 



-2 

^Aci (0= JP E 7C(0, T)bki e^c+jto,){iT+x,)^^{c+jt^,)T _ j ^ (3 

1+7-r 



(3-6) 



k=-MI2 l=^L 

kfH 



M/2 



k=-MH ^+ 7 — 

/fcjtO 



M/2 



+ r7P E 7t(Tt,7)6t-i (3- 

L . - - 1+y-^ 



/t=-M/2 

fctO 




+ 7? E 

k^MH ^+7“ 

+z?+zf,(=zf,,_.+z^;,,) 



zj — > 


(3-1) 


z?^ 


(3-2) 




(3-3) 




(3-5) 



(3- 
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IK WORST-CASE ANALYSIS 



A. Worst-case for S(t). 

T 

B. Worst-case for £ \S(t) dt - output of integrator. 

C. Special worst-case analysis. 

D. Special worst-case at the output of the integrator. 

A. Worst-case for S(t) 

A-1 Worst-case, ISI : = A_, L^ = oo 

From (2-2), we have : 

SJsM = jp b- - 1) X e<‘^T 

h=\ n - -I 

I I < 1 

= JPb-e-^‘ 0</<r (4.,) 
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A-2 Worst-case ACI : j = = b, Z, = ©o 



^ MU -2 

S7cA» = JPb I, {ic(0, 7)^7Tf - 1 )[ 1 



k=-M/2 

Jb* 



1+;^ 



/=- 



+7? 6 1; {JC(0, Tt)-is- } 

k i+y-r 



+7? 6 Z T)-^ _ e-^(«») g-<«yto*)7-) 

i i+y-r 



+ jc(x*, 7)-L- (eJ“‘k('-w _ e-c(i-x,)y } 

1+/-T 



= 7p 6 E IC(0, 7)- ^ g ( Ar) ^~{ j 

k 1+7^ 

+ Vp ^ X 71(0, X;t)-W g-(^70),)7’) 

;t 1+7^ 

+VP Z? X 7t(T;t, 7)-^ (g7to.(/-x.) _ g-c(/-x,) ^-(c+/(0,)7’) 

/t 1+7-# 

= JPbJ^ 71(0, 7)-^ g-<c>t/co*)r) 

/t 1+7^ 

+^^Pb X 71(0, 7)— T- (e7^('-^^> - 
k l+7-r 



= 7?i 1 7X(0, 7)^ 

it 1+7T 
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T 



B. 



Worst-case for J dt - output of inte 2 rator 



Output of photodetector : £ \ S (t) 

\sm= \SJt)+SJt)£S,,,(t)\^ 

= \SM + \SisM + (t) \' + (t) S]si(t) } 

+ 2 Re { Sg(t) S^(^j(t) } + 2Re { Sj^j(t) S^^^/t) } (4-3) 

B-1 From (2-1) 

(t) r = \SM = PKo + e-^V, 0<t<T, (4-4) 



r 



J \SB(t)Pdt = PTblo[l-^(l-e-^^) + ^0 



(4-S) 



0 



B-2 From (4-1) 



\S7si(t) r = Pbf 



0<t<T, (4-6) 



T 




(4-7) 



0 
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B-3 From (4-2) 



M/2 M/2 

= p\b\^ E E 

k=-M/2 m=-M/2 

h^O m^O 



^oif^(t-Xk) ^-ya)/n(r-x/n) 

(i+;?)(i-m 



= Pl6PEE 

k m 



(i+;^)(i-/=^) 



0<t<T, 



f 



l5;j,(0Pdr =Pl6|2 



ss 



(i+y?)a-/=^) 



T 

0 




cr 



sz 






jt W y(c cXl't/c)(^7c/ 



(£;<<• 



k ^m)T 



1), 



P7l6P E 



l+(-^)^ COi — 0/n 
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B-4 



2Re {SbH) SJS; (/) } = 2Fio.o6-(e-^' -e'^') 0<t<T 

j 2Re {Sb(0 SJ^; (0) dt = ^^^(1 - e~^^f 



B-5 



2Re {Sg(t) = 2Pbo,oRe {b (1 - e-")Z > 



= 2Pbo f)Re {b E ^ 0<t<T 

k i+^-r 



\2Re{Slit)S7amdt = 2^^. E ^ (^- ^)> 

A k C J C ^ J C 



B-6 



2Re S7cm = 2Pb.Re {Z. E ^ o<t<T 

k 



\2Re {S-;^{0S:tai‘)}dt 
0 



2PTb 

cT 



-Reibl, 

k 



e~J^k^k 



(1 - 



(4-11) 

(4-12) 

(4-13) 

(4-14) 

(4-15) 

(4-16) 
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C. Special worst-case analysis 



coyt = 2%kI/T, I : positive integer 

Tk = 0 

bo, = b-e{0,l} 
bkj = bk,o = b e {0,eJ^} 

Re{b}e {0, cos(|) } 

I6P6{0,1} i 

T 

Cd J \SB(l)l^dt = PTblo + (4-5; 

0 

Q2 J \SJi,(l)\ldt = ^ (1 - e-^^) (4-17 

0 



C-3 in (4-9) 

T 

J ^ X 0, (0A:?tC0;„ 



\\S7cMldt = PT\bV^^^ 



*^2 



k 



C-4 T 



J 2Re {Sb(()S';I; ( t) },,dt = (1 



QS j 2Re {5^(05:J^/(0 )scdt = (e-^r _ ^ 

n k ) 



C-6 



±1 2Re {Sj^;(t)S7a(t) }scdt = — 



(4-18) 



(4-12) 



(4-19) 



T (4-20) 
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D. Special worst-case at the output of the integrator 



1 

AT/ = J \S'’^(t)\ldt 



= ^ [ all terms in section C]. 



{ ii. ti - 



2 ( 



+ 



bl 
2cT 

bo.o b- 
cT 



cT 
,-2cT\ 



2cT 



(l_e-c7')2 



H\b\^ + ^^^Re{b}(b--bofi)-\'L-^,} 

k 



(4-21) 



46 



V. EXACT ANALYSTS - GENERAL CASE 



Now, we are to encounter the most complicated situation, for clarity, we 
first list the construction talbe, for which we are to follow : 



Sjsfi) S/t) 

Table 5-1 

I 

_| C A 

I I I 

2Re{S’,s,it)S^a(n} 2Re{S;(t)S^a(0} 2Re{SB{t)S}sm 

I 

I 

C-6 ^C-5 C-4 C-3 C-2 ^C-1 

I I I III 

2Re{Zl^Zr,} 2Re{ZlzZ} 2Re{Zlzh IzJI^ 

I 

I 

C-3-3 ^C-3-2 C-3- 1 

I I I 

misJzl 2Re{Zi_^zll,} IziJ^ Izlj^ 

We first calculate each term in Table (5-2), then integrate it. For each term 
belongs to we also discuss two cases: (0jt = 0)m and (0^;tc0m . 



^Ac/0 



yT yT 

^Zic ^0 ^0 



"c*,0 



7^ 



D 
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Since the computation is very complicated, here we just write down 
the results for each term, and omit the intermediate steps. 






(4-5) 



B. From (2-2) 



\S,si(t)\^ = 






j=—Lo 



0<t<T 



(5-1) 



-1 



C Calcualte \S - refer to (3-6) 



£=L 



IzJP = />e-2"jt(0, 7)1 E E bkj _ j)|2 (5.3^ 

k I 
T 

\\zl\^dt = ^J\- J e ic+JOik){lT-FZk) {c+j(i)ic)T — 1 ) 1 2 (5-4) 

i k I ' 1 +/-^ 
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CzZ Ts = min {Xk^m) 



IzJ' P = P S E lt(0, a; 

k m 



(i+;^)(i-7^) 



^g(c+;cOi)/ _ g(c+ya)i-)(-r+Xi )^ 






= pXXit(o,T.)6,,_,6;,., 



/W 






j-g;((Di_(D„)/ _^-(c+y(0„K g(c+y(0t)(-7’+-Ci) 



-{c-j(Sik)t Q^C-jdim ){-T-yZm ) (c+>C0fc)(-7’+Tt)g (c-/(Om )i.-T+T.m ) ] 



(5-5) 



I 



jlzJ‘lV/= JlzJ‘lv< 



(a) k^m, 



= -^ZZ^/r-l^w-l 






k m 



(l+7^)(l-;=) 









(^Q-{c+j(iirr,)ls _ 2^ g(t>»-ya>*)(-7’-Ki) 



1 , 



j L_ (^Q-iC-Mls _ 2^ Q{.C-j(H„){-T-^m) 



\-j^ 



.i _ 2) g(‘^7Wi)(-^-+^i)g(c-yoj„)(-7’+r„)j 



(5-6a) 
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(b) /: = m,C0A: = cOw 



f S ^ ICT, + - i)e(^->*X-7-«») 

Ic c ) I't/ c 

_l L_ (^Q-{c-j(^k)T^k _ \-^Q{c-j^k){-T+Xk) 

1-;? 



_l(g-2cXfc _ g(c+;COA:)(-7’+T)t) ^(c->*)(-7’+Tfc) 






cj- Y |4<?)2 



+2^e{ > (e -e(">*X- r+T.))}] 
l+7-r 



:5 5 /erm^/or IZj^l 

C-3-1 Xg = max(Tyt, Tw) 



m 









_g-(C+>fe)7’^ 






(<af) k^m, CDyt ^ CD 



W 



_ ^ y y /.-2cxg _ ^-2cr\ 

^ 7 ■ 



cT ^ ^ 2 

k m 



[g(c+;CO*:)TA:('l _ ^-(c+>;t)7’^ g(c-;<»n,)Tm Q - Q-{c-j<S> „ >")] 



(5-6b) 



(5-7) 



(5-8a) 



50 



(b) k = m, (Hk = dim 



k l+(-r)^ 



C-3-2 



3fc ^’’A^k‘^k~^rn‘^rn) 



\ZlJ^ = P S S 7l(Tj, 7) h,,bl^, 

Ic m U“t/ c )\^ J c ) 

|-g7((0 *-©„)/ _ Q-{c+j<a,n)t^ ^{c+j(s>k)-^k _ ^_(c_ytOfc)r ^(c-ycom)x. 

-\-e e ] 

T T 



\\zi;j^dt=\\zij^df. 



(a) k^m, (dk^a 

=§es6„6: 



/77 



PT , 3f( 

k m ■ (i+/^)(i-y^) 



cT ^ ^ .“m^ 



[ 



y(?-^) 



(^^y((o*-<Om)7’ _ ^y((Ofc-co;„)Tg^ 



^_£Z^ (^-(c+yo>;„)7 _ g-(c+yco„)Tg) 



1 +/ 



^^-ip-j(Sik)T _ g-(c-yco*)Tg^ 

1“^“ 

_i(g-2c7’ _ g-2cr^^ ^(c+y©*)T*g(c-> j 



(^) /: = w, coyt = CO 



/W 



= ^ S ^ .[cr- CT, + i(l - 

1 +;— 



(5-8b) 



(5-9) 



(5-lOa) 



(5-lOb) 
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C-3-3 



2Re{Zl_^Zl , } = IPReO, E nix ,,T)b 

k m 



^-jliO/.Xj.-tam'^m) 

(1+7%(1-^) 



^Q-{c+j(Sim)t Q{c+jVik)lk(^\ _g-(c+/®/t)7’^ 






A72 



* e 






= ^Re(L'Lbk-\b „. — — 

r 1 (p-ic+J(^m)T:g _p-(c+Jd)„)T\ 

Hi+jrp ^ 



_g-(c+y®jfc)7’^ pic+j(Ok)-Ck 



_l(^-2cTg _ ^-2cr^ j 



(5-11) 



(5-12) 
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C-4 



-2 



2Re{Zlzl‘ ]=2PRe>X'Z E 



k m l=-L 






^^-{c+j(n„)t^{c+j(iik){lT+Zk)^Q{c+j<i>k)T _ 



-let g (c-y© (C+/CO k)ilT+ik) (c+/(0 ;t)r _ j ^ j I 



\2Re{Z^Z]^]dt = J 2R6{ZqZq }dt for k?^:m and k=m 

0 0 



-2 






yt w l=-L 



(l+7^)(l-y^) 



rii-f iii p{c-J(i>m){-T-^-^m)’\\ 

L 2 ^ 
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C-5 



2Re{ZlzV,} = 



2PRe{Y, Z S 

lbk,ibl -i _e-(^->-)7') 

+6t,/6;o 



1 1 

j 2Re{ZlZl }dt = J 2Re{ZlzZ }dt 



-2 



= ^^e{SS E (e(cH.>.)7-_ i)g(cH->,)(/r-n>) 



yt /« /=-i (1+7-#)(1-;T) 






+ ^g-(c+7C0;„)T;„ -Q-ic+j^m)T>^ 






_ (g-2ctm _ g-2c7'^g(c-;(0„)T;„j j 



jy(c+;cOi.)(/r+i;) 



(5-15) 



(5-16) 
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C-6 



2Re{Z]^Zl^*} The product is only survived when Tm <t<Tk 



= 2PRe{Y.Y.Ti{'^m,Tk) 

k m 



(i+^fxwW' 



-bk-\b*m-\^ -2c/^ (c+y© feX-r+Tfc) ^ (c-y© m)Tm ^ ^ -(c-> 

+6t,.,6;o 

-bk-\bl,o 

-bk-\bl , 0 e-(">”)'e(">‘X-r+Ti,) 



(5-17) 
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f 2Re{Zl^Zll }dt = I 2Re{Z\‘ZV, }dt 

n 



2/^7’ D i V V 

cT {2^ 2^ n i,!^vi 

k m (^'t/ c )(1 y c ) 



r ^k,-\b^^ \ / 7 ( 

1-;^ 









7(“c c“ ) 

_ Q-{c-j<Sik)-Ck^Q(,C-j(Sirr,)-Cr 



i-y^ 



i^li£^(g-(c+;COm)Xm _ g-(c+7<Bm)XA:^g(c+yC0)t)(-r+TA:) 



1+y^ 



_g-2CT*:^^(c+;C0;t)(-7VTAr)^(c-yC0;„)T;„J J 



T T 

D. J 27?e{^s(/) 5;a(7)} = J 2^s(7) 5y5/(7) 

0 0 

-1 r 

= 2PZ)o,o E bo,i J dt 

i=-Lo 0 

i=-Lo 



(5-18) 



(5-19) 



I 
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E. 2Re{Sl(t)SACi(t)} 



= 2P/?e{ZZ 7t(0, 7)(e-" l)e(«-'“‘K'r+n) 

k I 1+y— 

+S it(0, u)bo 0*i -I - (£-" - e-2«)e(«>»X-r+tj)] 

A: ’ ’ I+y-r 

+1 7i(T*, 7)6o,o6*,-i^(e-"' - 

k 1+y— 

+Z n{Xk, T)bo,ohfi - e-(^->*)' - (g-^' - } (5-20) 

k ’ 1+y— 

T 

\lRe{Sl(t) SAcmdt 



_2 

'2J’T ry fV* V ^o,o^Jt,/(l-^"‘^0^ 

2 i+v^i 

K / C 



(g(c+y“2r)7’ _ g(c+ytat)(/7’+Tfc) 



+Z io.o6*,-, ^ [(^) - 

)t 1+2— j— i-y— 

-XS 6o,o6t-l^ ^^ 2e-k-2e-^r (j_ 

/t 1+2— 



+2 j OQ^obkO 



- IV? 



CD I, 
2 — 



/ “ V , n ) 

^ y c 



-( 



2 e"^*r- 2 e-^^-e- 2 ”fc+e- 2 ^ 3 ' 



g(c+ycofe)T*jj 



(5-21) 
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F. 2Re{S^siO)SAcm 



-1 



= 2P[ S 

i=—Lo 



Re{I. 1 7t(0, T)bu e-2" 

kl=-L ' l+;-r 



+E 7t(0, X*) bk,-\ ^ _ g-2«g(«>t)(-r+,t)) 

l“y c 



+S 7t(xi, 7) bt,-i ^ e-^'(l - e-«^x7-»)X)e(c+>t)x. 

k ^+J- 

+Z xt(xt, 7) bt,o ^ _ g-2c/g(c+>.)xc) } 

yt 



1 

\ 2Re {5;j/(0 SAci(t) } dt 



= ^ 2 bo,i - e'"^7?e{l2 



-1 



-2 



/=-Io 



yt / 2(i+y?) 



(^Qip+j<dik)T _ J^g(c+>*)/7’^ 



^ 1 “A: 

yt 1+7^ 



|-^ g *^^ >1 ^-{c+J<sik)T j 



1-7? 



+ E 



bk,-\ , 

” 2(1+;?) 



-CTa: _ ^CTfc^-2c7’^^J _ ^-(c+>/t)7’^ 



+ E -% [— W [e'"‘ - e-7“*x,g-(<>->»)r) _ _ gCt^-icT-. j j 

yt 1+7— 1-7— ^ 



(5-22) 



CXk 



(5-23) 
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VL EXACT ANALYSIS - SPECIAL CASE 



Now, we pul two constraints into all terms of exact analysis to form a 
special case ; 

1 . CO jt = I : positive integer. 

2. Tic = T for all adjacent channel k. 

For the purpose of normalization for future use, we set a = j , or Xk = x = al 

0 IT 

so that cO/tiA: = cOjtT = -^ar= 27c/:/a 
Still, we follow Table 5-2. 

A. 

B. 

C. 



]\SB{i)\^dt = PTbi^[\ 



2(1 (1-e-^^^) 
cT 2cT 



(6-1) 



0 

T 



0 i=-Lo (6-2) 



Special case for ACI 



C-1 



1 



k i=-L 1+y^ 



£2^(,2acr_,2,a-l).r)(,cr_ 1)2 I if 



(6-3) 
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C-2 



T 



T 



j \r^\^dt = \\Zl\^dt: 



0 



0 



(a) m, Cdk ^ CD/« 

= {cT)^PTYdYd 



^ ^ {cT+jlTikI){cT-j2TimI) 

I _g~y27i/a(A:-/w) 

jlnlik-m) 

^ (a- 1 )cT-^jlnmloi 

cT+J2nmI 

^ (a-l)cT-J2Tikla 

cT-jlnkl 

^2(a-\)cT_^-2cT 



IcT 



] 



(6-4a) 



(b)'/: = w, coyt = co,« 



= {cT)^PT E 



{cT)^+Q.nkrf‘ 




A: 



+2/?e{ 



^ -c7_^ (a- 1 )cT+j2 nkl a 



}] (6-4b) 



cT+jlnkl 
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C-3 3 terms for IZ^P 



C-3-1 ]\zf\^dl = \\zl_^\^dt 



(a) k ^ m,G)ic ^ 



m 



= \cTPT(\-e-^)W 2 



(b) k = m,(dk = 



m 






IV-il 



^ (cT)^+(2nkI)^ 



C-3-2 



I \Z^Q\^dt = ^ \Z^^^\^dt 



(a) m,G)k ^ (a 



m 



= {cT)^pt'Z'L 



bkfib*„fi 



k m 



{cT+j2nkI){cT-j2‘KmI) 



■{e 



-jliiJaQc-m) __ 



1) 



jlnlik-m) 



+ 



+ 



(a~ 1 )c 7 ’+/‘ 27 t/w/a _ 1 ^ 

cT+jl%ml 

{a-l)cT-j2nkla _ | ^ 

cT-j2%kl 



^(I_e 2 (a-l)c 7 ’) ^ 



2cT 



(6-5a) 



(6-5b) 



(6-6a) 
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(b) k-m, CO A: = CO 



m 



= {cTfPTY. , [1 - a + '' 



“ (cTf*Q.%kI)‘ 



IcT 



.2Re ] 



C-3-3 



llRe {Zl.iZl^o }dt = l2Re {Z^.^Z^}, }dt 



(^) m, (£ik ^ 



© 



m 



= 2(cT)^PTRe {HE 



k m 



6fc,-i6;,o (1-e-^^) 

{cT+j2nkr)(cT-j2Timr) 



^(^l_g(a-l)cT+j2nmJa^ ^ ^2(a-l)cT_i 

L cT+j2nml ' ^ 



2cT 



(b) k=m, ©A: = © 



m 



= 2{cT)^PTRe {E 






A- (cT)^H2nkI)‘^ 

|.(l_e(a-l)cr+y2it«a) ^ ^2(a-l)cr_j 

*- cT+j2itkI ' ^ 



2cT 



(6-6b) 



(6-7a) 



(6-7b) 
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C-4 



T 

J 2 Re { Zjz;* }dt = 



0 



\lRe { ZjZj’ )dt 



0 



(^) m, &k dim 



-2 

= 2 {cT)^PTRe {EZE 



k m I 



{cT+jlnkI){cT-jlnmI) 



^acT+j2nmJa_Y 

cT+jlnml 



+ 



g-cT_^(2a~\)cT 

2^ 



]} 



(b) k = m, ©A: = CO 



m 



= 2(cT)^PT {S S 



l 



k I 

gOjcT+J2iMa_^ g-cT_g(2a-l)cT 



cT+j2nkl 



2cT 



]} 



C-5 



1 1 
j 2 Re { Zjzf }dt J 2 Re { Zjzf }dt 



-2 



= 2 (cT)^PTRe {ESE 



(M)cT_^lcT 



ffi I (.cT+j2nkI)(cT-j2nmI) 2cT 

bubU iw6;.o 



+ 



cT+j^Titnl 



2cT 



]} 



C-6 



1 T 

{Z 5 Zr}c//= \ 2 Re {Zlzr}dt =0 



(6-8a) 



(6-8b) 






(6-9) 



( 6 - 10 ) 
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(6-11) 



D. ]2Re{SB{t)S]s,{t)}dt 



PTbo,o 

cT 



-1 



2 ^0,1 

/=— 2o 



1 

E. j 2Re{ Seif) 5^c/(0 }dt 



= 2/>r { S E 

k I 



2(cT+j2iikI) 



F, 



+2 

k 

+2 

k 



bo,obk-\cT (i_e-«cr) 2 g(a-i)cr 



cT-j2nkI 



cT+j2itkl f y27r^7 

2{cT+j2nkI) 



2cT 



] 



V bofibk,ocT |.(e-^2;tWa_j) ^^-cr-y2jcWa_g-acr) ^ 

cT+j2itkI I- ;27t/t7 ?7r J I 



cT-j2nkl 



I 

j 2^e{ 5;s,(0 SAcif) }dt 



-1 



= 2PT[ E 

/=— 2o 

{2. 2 

it / 



+E 

At 



2(cT+j2nkI) 

y-j2iikIa_^-ojc 

^ cT+j2nkI *■ cT-j2Ttkl 
bk,-\ (1-e-^^) 



6*-lCr ^(g-/2JtWa_g-acr) (^(a-\)cT _^-(a+\)cT>^ ^ 

' ^ L ^T-iOtrEJ 2 cT -I 



2(cT+j2nkI) 



V (g^r_^-cr-;2nA:/a) (g-acr_g(a-2)cr) ^ 

+Zi L cT-j2nkl ~ '>rT J / 



cT+j2%kl 



2cT 



2cT 



(6-12) 



( 6 - 13 ) 
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APPENDIX B 



Programs 



% exact analysis with optimal threshold signal out of the integrator, the signal X contains desired 
signal and ACI & ISI and postdetect noise ,The formula for the BER is : 



Where = 



2% T 



[J-J I JQ( 

MO MO 






)d^ -MU . . . d^M/2 dx -MU ...dxMH 



1% T 



+ J...J J...J ^)d^-MI2^’-d^Mll dX-MIl—dXMll] 

AT 0 AT 0 



M is the number of adjacent channels. Here we assume all channels have the same phase (j) and 

time delay T i.e, (j)jt = Tjt = T hence for our model the power of (1/271) and (1/T) are fixed 
to 1, therefore we only have two integrals and two arguments. 

M=4; k=[-M/2:-l l:M/2]; 
m=k; 

% produce the controlled matrix b to control 64 different bit patterns 
ml=[ zeros(l,32) ones(l,32) ]; 
m2=[ zeros(l,16) ones(l,16) zeros(l,16) ones(l,16)]; 
m3=[]; 
for i=l:4 

m3=[m3 [zeros(l,8) ones(l,8)]]; 
end 
m4=[]; 
for i=l;8 

m4=[m4 [zeros(l,4) ones(l,4)]]; 
end 
m5=[]; 
for i=l:16 

m5=[m5 [zeros(l,2) ones(l,2)]]; 
end 
m6=[]; 
for i=l;32 

m6=[m6 zeros(l,l) ones(l,l)]; 
end 
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b=[m 1 ;m2;m3 ;m4;m5;m6]'; 



% signal to noise ratio range in dB 
RPTN_DB=[10:0.5;25]; %RPTN_dB 
ppp=l 0.^(0. 1 *RPTN_DB); 
lenl=length(ppp); 

% solalph function is provided by Professor Randy L. Borchardt 



n=10; 

[bpx,wfx]=grule(n); %bpx=bpy ,wfx=wfy 

% single channel 
BER0=0.5*erfc(ppp/8'^0.5); 

pp=[]; thresh=[]; %thresh is not normalized threshold 

for CT=[5 10 15 20] % cT is Fabry-Perot filter parameter 

ifCT=5 

qqq=[2;12 17 20]; 
elseif CT=10 
qqq=[3;12 17 20]; 
elseif CT= 15 
qqq=[4;12 17 20]; 
elseif CT=20 
qqq=[5:12 17 20]; 
end 

pe=[]; threshl=[]; 
for I=qqq 

BER=[]; thresh2=[]; 
for RPTN=ppp 

ap=linspace(0, 1 82, 1 1 ); % approximated thresholds 

x3=0; % first few loops to find out the threshold which make the 

% BER minimun don't care about the scale 

for i=l :64 x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 
solalph('xxl r,0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 

lef=ap(ind)-16.6; 
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if lef<0 ,lef=0; end % to aviod the threshold go beyond the negative side 



ap=linspace(lef,ap(ind)+ 1 6.6, 1 1 ); 
% two more time to find alpha 
x3=0; 



fori=l;64 x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 
solalph('xx 1 1 ',0,2*pi,2,bpx,wfx,0, 1 ,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[ val, ind]=min(x3 ); 
lef=ap(ind)-3.1; 

if lef<0 ,lef=0; end 

ap=linspace(Ief,ap(ind)+3 .1,11); 

% three more time to find alpha 
x3=0; 

for i=l;64 
x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 
solalph('xxl r,0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 

lef=ap(ind)-0.57; 

if lef<0 ,lef=0; end 

ap=linspace(lef,ap(ind)+0.57,8); 

% four more time to find alpha 
x3=0; 

for i=l:64 
x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 

solalph('xxir,0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 

end 

[val,ind]=min(x3); 

ap=ap(ind); 
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% after find optimal alpha use double integration 
qq3=0; 



for i=l :64 

qq3=qq3+( dbgquadm('xx00',0,2*pi,2,bpx, wfx,0, 1,2,.. 
bpx,v^,ap,b,CT,i,I,k,m,RPTN)+... 
dbgquadm('xxir,0,2*pi,2,bpx,wfx,0,l,2,... 
bpx,wfx,ap,b,CT,i,I,k,m,RPTN) )/(256*pi); 
end 

BER=[BER qq3]; thresh2=[thresh2 ap]; 

if qq3<10'^(-15) %for save time only interesting in 10^(-15) 
ii=find(ppp=RPTN); 

BER(ii+l:lenl)=5*10^(-l 16)*ones( l,length(ii+l:lenl) ); 
thresh2(ii+l :lenl)=(ap+2)*ones(l,length(ii+l :lenl) ); 
break 
end 



pe=[pe;BER]; thresh l=[threshl ;thresh2]; 



pp=[pp;pe]; thresh=[thresh;threshl ]; 
end 

time2=toc; 
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% worse case form appendix setting optimal threshold equal (x0+xl)/2 
M=4; 

I=linspace(0,6,101); % set I value in linspceinteger to find minimum I 

RPTN_DB=0:0.2;20; % x-axis dB range 

RPTN=1 0.^(0. 1 *RPTN_DB); % change to ratio 

% single channel 

BER0=0.5*erfc(RPTN/8^0.5); 

%find optimal alpha 
pp=[]; 

forCT=[5 10 15 20] 
x3=0; 

fork=l;0.5*M 

x3=x3+2./(l+(2*pi*k*I/CT).^2); 

end 

x0=(l-exp(-2*CT))/(2*CT)+x3*(l+2*(l-exp(-CT))/CT); 
xl=l -2*(1 -exp(-CT))/CT+(l -exp(-2*CT))/(2*CT); 
[val,ind]=min(abs(xO-xl )); 

I=ceil(I(ind)); 

I=(I(ind)); % our minimum I value 

% set different I and to find the BER 
ss=I+l:I+10; 

BER=[J; 
for I=ss 
x3=0; 

fork=l:0.5*M 

x3=x3+2/(l+(2*pi*k*I/CT)^2); 

end 

x0=( 1 -exp(-2*CT))/(2*CT)+x3 *( 1 +2*( 1 -exp(-CT))/CT); 
xl=l-2*(l-exp(-CT))/CT+(l-exp(-2*CT))/(2*CT); 
BER=[BER;0.5*erfc( RPTN*(xl-x0)/8'^0.5 ) ]; 
end 

pp=[pp;BER]; 



semilogy(RPTN_DB,BERO.’-’,RPTN_DB,pp(;,;)) 
axis([10 19 10^(-15) 1]) 
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